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ABSTRACT 


The  "Computer  Assisted  Target  Analysis"  program  cataloged 
in  computer  form  the  ground  truth  data  collected  at  the  UNDERBRUSH 
Test  Range,  Eglin  Air  Force  Base,  Florida.  A  Fortran  IV  program 
compatible  with  the  GE-645  computer  was  developed  for  storing, 
retrieving  and  correlating  the  data,  and  providing  a  means  for 
updating  as  more  data  are  compiled.  Contained  within  is  a  complete 
program  documentation  including  step  by  step  operating  procedures 
and  logic  flow  diagrams.  In  addition,  the  data  reduction  method¬ 
ology  and  program  test  results  are  described. 
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1.  INTRODUCTION 

The  "Computer  Assisted  Target  Analysis"  study  is  designed 
to  catalog  in  computer  form  the  ground  truth  data  collected  at  the 
UNDERBRUSH  Test  Range,  Eglin  Air  Force  Base,  Florida,  to  store, 
retrieve,  and  correlate  the  data,  and  to  provide  a  means  for  fu¬ 
ture  updating  as  more  data  are  compiled. 

Basically  this  report  can  be  divided  into  two  segments. 
The  first  part  presents  the  complete  program  documentation  for  the 
computer  program  and  associated  subroutines  which  were  developed 
in  support  of  this  study.  Each  salient  portion  of  the  program  is 
described  in  dc’ail,  including  flow  charts,  program  operation 
description,  inputs  and  outputs,  and  operating  instructions.  This 
documentation  section  is  intended  primarily  for  use  by  persons 
familiar  with  the  GE-645  computer  and  Fortran  IV  language.  The 
information  contained  in  this  documentation  will  permit  any  modi¬ 
fication  to  the  program  to  be  undertaken  in  an  efficient  manner. 

The  second  part  discusses  the  techniques  used  in  data 
reduction,  presents  the  computer  test  results  undertaken  to  simu¬ 
late  operational  requirements,  and  considers  potential  applications 
of  the  program.  Information  on  the  alphabetic  and  numeric  codes 
used  in  the  data  reduction  are  explained  in  Appendix  A. 
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2.  GENERAL  SYSTEM  DESIGN  AND  OPERATION 

The  program  created  for  the  "Computer  Assisted  Target 
Analysis"  study  was  designed  for  the  GE-64S  computer  using  the 
GECOS  III  configuration.  Fortran  IV  language  was  used  throughout 
the  entire  programming  cycle. 

The  basic  program,  mnemonic  COMATAN,  is  composed  of  the 
following  eight  subroutines:  STORE,  RETRIV,  BRITE,  CORREL,  ARRAY, 
MISR,  PLOT,  and  SCAL.  A  detailed  description  of  these  subroutines 
is  presented  in  subsequent  sections. 

Before  examining  the  component  parts  of  the  program  an 
understanding  of  the  overall  system  design  is  important.  This 
section,  therefore,  describes  the  operation  of  the  COMATAN  pro¬ 
gram,  the  options  available,  and  the  set-up  of  the  card  decks. 

In  order  to  operate  this  program  two  magnetic  tape  units, 
a  card  reader,  and  a  printer  must  be  available.  Of  the  two  magne¬ 
tic  tapes  used  for  this  study,  one  contains  all  the  encoded  UNDER¬ 
BRUSH  ground  truth  data,  and  the  other  is  a  scratch  tape  used  for 
retrieval  compilation. 

Basically  the  program  is  written  in  three  parts:  stor¬ 
age,  retrieval,  and  correlation  phases. 

2.1  Storage  Phase 

The  function  of  the  storage  phase  is  to  compile  all  the 
formatted  ground  truth  data  onto  a  master  magnetic  tape.  There 
are  three  options  available: 


1.  Initial  Store 

2.  Additional  Store 

3.  Housekeep  Store. 

The  initial  store  option  pertains  to  the  loading  of  data 
from  punch  cards  onto  the  magnetic  tape  for  the  first  time.  This 
option  is  exercised  when  creating  a  master  tape  from  the  data  cards. 

Since  ground  truth  data  are  being  continually  generated, 
it  would  be  advisable  to  update  the  master  tape  periodically.  The 
second  option,  additional  store,  provides  for  this  operation.  Here 
additional  data  in  the  form  of  punch  cards  can  be  added  to  the 
master  tape. 

The  last  option,  housekeep  store,  is  used  for  a  tape-to- 
tape  transfer  of  data  from  the  master  tape  which  may  become  worn 
to  a  fresh  tape.  It  serves  as  a  means  of  preventing  the  loss  of 
data  on  an  old  tape  and  having  to  generate  a  new  tape  from  the 
punch  cards  again. 

2 . 2  Retrieval  Phase 

The  retrieval  phase,  wherein  all  the  data  are  searched 
and  the  desired  information  printed  out,  is  the  crux  of  the  entire 
program.  Six  options,  two  input  and  four  output,  are  available. 

2.2.1  Input  Options 

1.  The  requestor  may  rup  up  to  10  cases  with  each 
case  composed  of  up  to  35  limits. 
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2.  The  above  option  may  be  run  for  any  number  of 
users  on  the  same  execution.  This  allows  for 
an  almost  indefinite  number  of  parameter  vari¬ 


ations  to  be  compiled  from  the  master  tape. 

2.2.2  Output  Functions 

1.  Retrieve  the  requested  records  on  input-for¬ 
matted  punch  cards. 

2.  Retrxeve  the  records  with  appropriate  case 
number  on  a  blocked,  input-formatted  magnetic 
tape . 

3.  Retrieve  the  records  with  the  appropriate  case 
number  to  an  on-line  printer. 

The  four  output  options  are:  neither  1  nor  2,  1  not  2,  2  not  1, 
and  both  1  and  2.  In  all  instances  the  retrieval  records  are 
printed  out,  function  3. 

2 . 3  Correlation  Phase 

There  are  no  options  in  the  correlation  phase,  instead 
the  retrieved  data  are  processed  through  the  entire  correlation 

program. 


2.3.x  Correlation  Phase  Input 

Input  to  the  correlation  phase  is  from  the  magnetic  tape 
generated  during  the  retrieval  phase,  limited  to  only  one  user, 
but  with  up  to  ten  cases  available.  A  control  card  stating  each 


individual  case  number  to  be  used  and  the  respective  X  and  Y  pa¬ 
rameters  to  be  correlated  must  be  provided. 

2.3.2  Correlation  Phase  Output 

Output  to  the  correlation  phase  consists  of  the  follow¬ 
ing  items: 


1. 

X-mean 

2. 

Y -mean 

3. 

X-variance 

4. 

Y- variance 

5. 

Correlation 

coefficient 

6. 

Regression 

line  plotted. 

2.4  Card  Deck  Set-Up 

The  COMATAN  program  has  five  execution  options  which 
are  presented  below. 

1.  Store  "N 

No  Retrieve  ^1,  0,  0 
No  Correlate 


2.  Store 
Retrieve 
No  Correlate 


>1,  1,  0 


3.  Store 
Retrieve 
Correlate 


1,  1,  1 
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!'  V 

»T 


4. 


No  Store 
Retrieve 
No  CorrelateJ 


1,  0 


5.  No  Store 
Retrieve 
Correlate 


jo.  1,  1 


On  the  following  five  pages  each  of  these  execution  options  are 
illustrated  with  respect  to  the  card  deck  set-up. 


2. 5  Data  Cards  and  Control  Cards 

In  order  to  utilize  the  COMATAN  program,  as  in  the  case 
with  any  computer  program,  not  only  must  the  card  deck  be  properly 
sequenced,  but  the  various  cards  must  be  correctly  keypunched. 

This  subsection  tabularizes  the  control  and  data  cards,  showing 
the  required  format,  and  describing  various  encoded  parameters. 


2.5.1  Execute  Option  Card 


COLUMN 

1 

2 

3 

FOR  EACH 


ITEM 


FORMAT 


Storage  Phase  Flag  (ISTORE)  II 
Retrieval  Phase  Flag  (IRETRV)  II 
Correlation  Phase  Flag  (ICOREL)  II 
J"  0  ■  Do  not  execute  phase 
1  -  Execute  phase 
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CARD  DECK  SET-UP:  EXECUTION  OPTION  1 
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CARO  DECK  SP.T-UP:  EXECUTION  OPTION  2 


CARD  DECK  SET-UP:  EXECUTION  OPTION  3 


2.S.2 


Storage  Option 

Card 

COLUMN 

ITEM 

FORMAT 

1 

Type  of  Storage  to  be  done 
(IRUN) 

11 

0  •  Initial  data  store 

1  ■  Copy  tape  and  add  data 
from  cards 

2  ■  Copy  from  old  tape  all 
data  onto  new  tape- 
"Housekeeping" 

Blocking  Indices  Card 

COLUMNS 

ITEM 

FORMAT 

1 

Tape  Drive  Number  Where  Data 
Tape  is  Located( ITAPE) 

11 

2  -  6 

Number  of  Full  Blocks  of  Data 
On  Tape(IPABCl) 

IS 

7  -  8 

Number  of  Records  in  Last 
Block  ( IPABC2) 

12 

NOTE:  A  Full  Block  of  Data  ■  Ten  Records 


A  Record  ■  Two-Card  "UNDERBRUSH"  Data  Sets 


2.5.4 


Data  Deck 

CARD  1 


COLUMNS 


1 

-  5 

MISNAM 

6 

-  9 

MISNUM 

11 

SENTYP 

13 

-  14 

SENDES 

16 

-  21 

I  DATE 

23 

-  25 

AIRCRF 

12 


ITEM 

FORMAT 

Mission  Name 

A5 

Mission  Number 

A4 

Sensor  Type 

11 

Sensor  Designation 

12 

Date 

16 

Aircraft 

A3 

27 

LOCAT 

*  Location 

A1 

29 

-  32 

SUNRIZ 

-  Sunrise 

14 

34 

-  37 

SUNSET 

-  Sunset 

14 

39 

-  42 

MUNRIZ 

-  Moonrise 

14 

44 

-  47 

MUNSET 

-  Moonset 

14 

49 

-  52 

LOCUM 

-  Tine  (Local) 

14 

54 

SKYCOV 

-  Sky  Cover 

11 

56 

-  57 

CLHITE 

-  Height  (Clouds) 

12 

59 

CLTYPE 

-  Type  (Clouds) 

11 

61 

-  62 

WINDIR 

-  Wind  Direction 

A2 

64 

-  65 

WINSPD 

•  Wind  Speed  (Knots) 

A2 

67 

-  68 

WESTAT 

-  Weather  Status 

A2 

70 

-  72 

WETBUB 

-  Wet  Bulb  (F) 

13 

**  77 

-  80 

Sequential  Card  Identification-  14 
Never  Read,  Written,  or  Processed 

CARD 

l 

COLUMNS 

ITEM 

FORMJ 

1 

-  3 

DRYBUB 

-  Dry  Bulb  (F) 

13 

5 

-  7 

RELHUM 

-  Relative  Humidity  t 

13 

9 

-  11 

SOLRAD 

-  Solar  Radiation 

13 

13 

PRECIP 

-  Precipitation 

A1 

15 

GRDVIS 

-  Ground  Visibility 

11 

17 

-  18 

DEWPTF 

-  Dew  Point  (F) 

12 

20 

-  22 

REFCUB 

-  Reference  Cube 

13 

24 

-  26 

BGRRES 

•  Background  Results 

13 

28 

-  30 

BGRTEM 

•  Background  Temp  (C) 

13 

13 


32 

MUNPHA 

-  Moon  Phase 

11 

34 

-  37 

SBRITW 

-  Standard  Brightness 
(Whole) 

14 

38 

-  42 

SBRITD 

-  Standard  Brightness 
(Decimal) 

15 

44 

-  47 

BBRITW 

-  Background  Brightness 
(Whole) 

14 

48 

-  52 

BBRITD 

-  Background  Brightness 
(Decimal) 

15 

54 

-  56 

SUMOSH 

-  Sun/Moon  Shadow 

13 

58 

FOLIAG 

-  Foliage 

11 

**  77 

O 

00 

• 

Sequential  Card  Identification-  14 
Never  Read,  Written,  or  Processed 

Cards  1  and  2  are  repeated  as  many  times  as  there  are 


records.  The  last  two  cards  of  the  deck  are  always: 

COLUMNS  ITEM  FORMAT 

1  -  5  Termination  Card  (Nines)  99999  1 

and  a  blank  card 

2.5.5  Retrieval  Option  Card 

COLUMNS  ITEM  FORMAT 

1  Retrieved  Records  on  Input-  II 

Formatted  Punched  Cards 
(IPUNCH) . 

2  Retrieved  Records  with  Case  II 


Number  on  Blocked,  Input- 
Formatted  Magnetic  Tape  (IKEEP). 
Required  for  correlation  input. 

J  0  •  Do  not  perform  the  option 
FOR  EACH  < 

1  1  -  Perform  the  option 


In  each  case  the  case  number  is  printed  followed  by  a 
one-line  print  out  of  the  retrieved  record. 
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2.5.6  Search  Parameter  and  Interval  Card  Set 

COLUMNS  ITEM  FORMAT 

1-6  Search  Parameter  (INAM)  A6 

7-9  Last  Card  Set  (IEND) -Right 

justified.  Used  only  after 
last  search  set  for  each  user 
group . 

12-21  1  Lower  Limit  of  Interval- 

Right  Justified 

12-21  2  Lower  Limit  of  Whole  Part  of 

Word-Right  Justified 

24-33  1  Upper  Limit  of  Interval- 

Right  Justified 

24-33  2  Lower  Limit  of  Decimal  Part 

of  Word 

36-45  2  Upper  Limit  of  Whole  Part  of 

Word 

48-57  2  Upper  Limit  of  Decimal  Part 

of  Word 

The  2  indicates  the  item  description  when  the  intervals 
are  represented  by  words  of  7  character  length  or  more;  such  as, 
standard  brightness  and  background  brightness,  each  of  which  con¬ 
tains  9  characters  (four  whole  number  digits  and  5  decimal  digits) 

The  above  card  represents  "one"  specification  of  a  pos¬ 
sible  35  available  for  retrieval  of  each  "Project  UNDERBRUSH" 
record  stored  on  magnetic  tape.  For  example,  a  typical  search 
may  involve  location  of  all  entries  in  June  1968,  between  1830 
and  2230  (local  time)  that  only  include  clouds  over  3,000  feet 
and  a  south  wind.  The  series  of  cards  to  describe  such  a  case 
are : 

15 


< 


WINDIR 

S 

S 

CLHITE 

3 

99 

LOCTIM 

1830 

2230 

I  DATE 

680601 

680630 

As  illustrated  in  the  card  deck  set-up  in  the  event  that 
multiple  retrieval  operations  are  requested,  a  blank  card  is  used 
to  separate  each  set  of  search  cards.  There  may  be  as  many  as  ten 
»uch  sets  within  each  user  group,  however,  the  last  set  must  be 
terminated  not  with  a  blank  card,  but  by  a  card  containing  char¬ 
acters  as  follows: 

COLUMNS  ITEM  FORMAT 

1  -  6  Terminal  Nines,  999999  A6 

9 


Last  User  Group  Flag  (IEND) -  13 

Right  Justified 


IEND  ■  0,  There  are  more  user  groups  of  search 
criteria  to  examine 

IEND  *  1,  There  are  no  more  user  groups  to 
examine,  return. 


CORRELATION  REQUEST  CARD 
COLUMNS 


ITEM 


1  -  2 

7  -  12 
14  -  19 
25 


Case  Number  of  Retrieved 
Records  to  be  used  in 
Correlation  (ICASE) 


FORMAT 

12 


Abscissa  of  Correlation  (IABS)  A6 

Ordinate  of  Correlation  (IORD)  A6 

0,  There  are  more  cases  to  be 
correlated-continue  in  the 
correlation  cycle. 
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1,  There  are  no  additional  cases 
to  be  correlated,  end  the 
computation. 


i 


3. 


COMATAN 


3 . 1  Summary 

This  routine  controls  the  main  logic  of  the  Computer 
Assisted  Target  Analysis  Program. 

3. 2  Restrict  ions /As  sumptions 

As  many  as  two  magnetic  tape  units,  a  card  reader,  and 

a  printer  must  be  available  for  use  of  this  routine. 

3 . 3  Options 

There  are  five  execution  options: 

1.  Store  only 

2.  Store  £  retrieve 

3.  Store,  retrieve  and  correlate 

4.  Retrieve  only 

5.  Retrieve  and  correlate. 

3.4  Method 

This  program  first  reads  a  control  card  of  three  flags 

to  determine  the  program  phases  to  use  on  a  particular  execution. 

Only  a  one  (1)  in  the  proper  place  designates  use  of  a  phase.  The 

first  test  is  for  use  of  the  storage  phase  of  the  program.  If 

affirmed,  a  message,  "MOUNT  USED  TAPE  ON  UNIT  1,  NEW  ON  UNIT  2," 

to  prepare  for  the  storage  routine  options  is  printed. 

NOTE:  "USED"  refers  to 

1.  A  blank  tape  onto  which  card  images  will 
be  stored. 
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2.  A  tape  containing  soat  card  iaagea  to  be 
copied  onto  unit  2  and  added  there. 

3.  A  coapleted  tape  of  card  iaagea  to  be  copied 
onto  unit  2  for  "housekeeping"  purpose!. 

Aftor  subroutine  STORE  is  called  a  Message,  "FINISHED  THE  STORAGE 

PHASE  -  DATA  ON  UNIT  _ ,"  to  identify  the  location  of  the  aost 

recently  stored  data  is  printed.  Then  the  retrieval  phase  is 
tested  for  use.  If  the  test  is  negative,  execution  stops.  If 
positive,  the  correlation  phase  is  executed  and  "FINISHED  CORRELA¬ 
TION  PHASE"  is  printed.  This  is  the  prograa  termination.  However, 
in  event  that  the  storage  phase  is  not  elected  for  execution,  a 
message,  "MOUNT  DATA  TAPE  ON  UNIT  1,"  is  written,  the  value  one  (1) 
is  assigned  to  the  tape  unit  naae,  and  the  retrieval  phase  flag 
is  tested  for  use.  Here,  if  RETRIV  is  not  elected,  execution 
stops.  If  it  is  elected,  the  argument  control  card  is  read  as 
input  to  RETRIV.  After  RETRIV  is  called  and  executed  "FINISHED 
RETRIEVAL  PHASE"  is  written  and  the  correlation  phase  is  tested 
for  use.  If  negative,  terminate  execution.  If  positive,  execute 
CORREL,  write  termination  message  and  end. 

3. 5  Calling  Sequence 
None . 

3 . 6  Common  Storage 
None . 

3 . 7  Subroutines  Required 
STORE 
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RLTRIV 
COR  Rill. 


3 . 8  Accuracy 

Single  precision  floating  point  word. 

3 . 9  Crror  Messages 
None . 


3.10 

fonaat  : 


Input  Format 

The  subroutine  option  card  must  be  rzad  in  the  following 


COLUMN  I THM  FORMAT 

1  Flag  for  use  of  storage  phase  II 

2  Flag  for  use  of  retrieval  phase  II 

3  Flag  for  use  of  correlation  II 

phase 


A  one  (1)  represents  use  of  the  phase;  a  zero  (0)  represents  non 
use  of  the  phase. 


If  the  storage  phase  is  not  used  and  the  retrieval  phase 
is  required,  a  card  must  be  read  to  define  the  calling  sequence  of 
subroutine  RLTRIV.  All  items  arc  to  be  right  justified. 

COLUMN  I THM  FORMAT 

1  Tape  unit  number  on  which  the  II 

"Project  UNDHRBRUSH"  data  tape 
(ITAPF.)  is  to  be  placed  for  use. 

Th is  i s  un i t  1 . 
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2  -  5 


The  nuaber  of  coapletely  filled 
block*  of  data  ( I P ABC  1 )  on  that 
tape.  14 

6  •  7  The  nuaber  of  two-card  data  sets 

in  the  unfilled  block  ( I PABC2 ) .  12 

The  values  of  IPAHC1  and  IPABC‘2  are  printed  during  any  execution 
of  (he  storage  phase  of  the  prograa.  These  values  aay  be  copied 
and  applied  to  the  above -aent ioned  calling  sequence  control  card. 

5.11  Output  Foraat 

The  following  on-line  printouts  may  occur  depending 
upon  the  mode  of  operation. 

1.  Mount  data  tape  on  unit  1 

2.  Mount  used  tape  on  unit  1,  new  on  unit  2 

5.  Finished  the  storage  phase  -  data  on  unit  _ 

4.  Finished  retrieval  phase 

5.  Finished  correlation  phase. 
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4.  STORE 

4 . 1  Summary 

This  subroutine  stores  "Project  UNDERBRUSH"  records  on 
tape  with  record  addition  and  protection  options. 


4 . 2  Restrict ions/Assumpt ions 

The  original  source  document  data  must  be  converted  into 
a  machine  readable  form.  The  encoded  data  for  each  set  must  appear 
in  two-card  sets.  One  or  two  magnetic  tapes  are  required  for  this 
subroutine  depending  upon  the  option  used. 


4 . 3  Opt  ions 

It  makes  provisions  to  store  the  data  initially  onto 
tape  from  cards  or  copy  the  previously  stored  data  from  an  original 
tape  to  a  new  tape  and  add  any  new  data  to  the  resultant  tape. 

Since  tape  records  are  so  easily  damaged,  all  data  stored  on  a 
tape  should  be  copied  onto  fresh  tape  periodically. 

4 . 4  Method 

This  subroutine  reads  the  storage  option  card  and 
branches  to  execute.  The  data  are  handled  in  blocks  of  10  two-card 
sets  for  optimum  use  of  tape.  At  the  end  of  each  option  the  num¬ 
ber  of  filled  blocks  (IPABC1)  on  tape  and  the  number  of  two-card 
sets  in  the  unfilled  block  (IPABC2)  are  reported.  The  tapes  are 
then  rewound  and  the  execution  continues. 

Provision  is  made  in  this  subroutine  to  distinguish 
between  a  wind  speed  reading  of  zero  miles  per  hour  versus  a  case 
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where  no  wind  speed  reading  was  made  (a  blank).  A  two-character 
Hollerith  letter,  IB,  is  used  to  designate  a  blank  wind  speed 
field. 

4 . 5  Calling  Sequence 

This  subroutine  is  called  by: 

CALL  STORE  (ITAPE,  IPABC1,  IPABC2) 

Where , 

ITAPE  ■  the  tape  unit  on  which  the  data  tape 
finally  appears, 

IPABC1*  the  number  of  completely  filled  blocks 
of  data  on  tape, 

IPABC2-  the  number  of  two-card  data  sets  in  the 
unfilled  block. 

4  .  b  Common  Storage 

None . 

4 .  7  Subroutines  Required 

None . 

4 . 8  Accuracy 

Single  precision  floating  point  word. 

4 . 9  Error  Messages 

Only  one  error  message  may  appear  for  this  subroutine. 
In  the  event  that  a  character  other  than  0,  1,  or  2  appears  for 
the  storage  option  control  (IRUN),  this  message  will  result: 
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’’STORAGE  OPTION  CONTROL  (IRUN)  CONTAINS  AN 


ILLEGAL  INTEGER  CONSTANT” 
and  the  integer  constant  will  be  printed. 

4.10  Input  Format 

Input  is  provided  to  the  routine  through  cards  or  tape 
or  tape  and  cards.  The  card  input  forma*  is  as  follows: 

CARD  1  -  STORAGE  OPTION  CARD 
.COLUMN  ITEM  FORMAT 

1  Storage  option  (IRUN)  Integer 

IRUN  -  0,  store  data  initially  -  read  and  store  initial 
data  from  card  to  tape  with  a  terminal  entry  of 
nines. 

IRUN  ■  1,  copy  tape,  add  data  -  read  old  tape  data  onto 
a  new  tape,  sense  terminal  entry,  add  new  data 
from  card  to  new  tape,  and  apply  terminal  nines. 

IRUN  ■  2,  copy  tape  *•  copy  data  from  old  to  new  tape  as 
a  "housekeeping"  exercise  for  protection  of  tape- 
stored  data. 

CARD  SET  2  -  ENCODED  "UNDERBRUSH"  DATA 

CARD  1 

COLUMNS  ITEM 

1  -  5  Mission  Name 

6  -  9  Mission  Number 


FORMAT 

AS 

A4 
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11 

Sensor  Type 

11 

13  -  14 

Sensor  Designation 

12 

16  -  21 

Date 

16 

23  -  25 

Aircraft 

13 

27 

Location 

11 

29  -  32 

Sunrise 

14 

34  -  37 

Sunset 

14 

39  -  42 

Moonrise 

14 

44  -  47 

Moonset 

14 

49  -  52 

Time  (Local) 

14 

54 

Sky  Cover 

11 

56  -  57 

Height  (Clouds) 

12 

59 

Type  (Clouds) 

11 

61  -  62 

Wind  Direction 

A2 

64  -  65 

Wind  Speed  (Knots) 

A2 

67  -  68 

Weather  Status 

A2 

70  -  72 

Wet  Bulb  (F) 

13 

CARD  2 

COLUMNS 

ITEM 

FORMAT 

1  -  3 

Dry  Bulb  (F) 

13 

5  -  7 

Relative  Humidity  % 

13 

9  -  11 

Solar  Radiation 

13 

13 

Precipitation 

A1 

15 

Ground  Visibility 

11 

17  -  18 

Dew  Point  (F) 

12 

20  -  22 

Reference  Cube 

13 

27 


24  -  26 

Background  Results 

13 

28  -  30 

Background  Temperature  (C) 

13 

32 

Moon  Phase 

11 

34  -  42 

Standard  Brightness 

19 

44  -  52 

Background  Brightness 

19 

54  -  56 

Sun/Moon  Shadow 

13 

58 

Foliage 

11 

There  may  be  as 

many  of  the  two-card  sets  as  necessary 

to  complete  the  data  available.  This  set  of  cards  must 

nated  by  these  two  cards: 

be  termi 

COLUMN 

ITEM 

FORMAT 

CARD  1  1-5 

99999  Termination 

A5 

CARD  2  1-80 

Blank 

4.11  Output  Format 

Output  occurs 

on  magnetic  tape  in  card  image 

blocks  of 

ten  two-card  data  sets  per  filled  block.  The  number  of  completPly 
filled  blocks  of  data  stored  on  tnpe  (IPABC1)  and  the  number  of 
two-card  data  sets  in  the  unfilled  block  CIPABC2'  are  then  printed. 
A  terminal  entry  of  9999  appears  at  the  end  of  the  records. 
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s. 


RHTRIV 


5 . 1  Summary 

This  subroutine  performs  a  search  of  the  data  stored  on 
magnetic  tape  to  permit  access  to  all  of  the  information  contain¬ 
ing  u  given  data  attribute. 

5 . 2  Rest  rict ions /As  sump t ions 

This  subroutine  assumes  input  from  magnetic  tape  of 
"UNDHRBRUSH"  encoded  data  records  and  presence  of  a  magnetic 
scratch  tape. 

5 . 3  Options 

This  subroutine  has  the  option  of  punching  all  retrieved 
records  for  one  or  more  users  on  cards  in  input  format  or  writing 
on  scratch  tape  all  retrieved  records  with  case  number  for  cor¬ 
relation  input.  A  documented  list  of  all  retrieved  items  will 
appear. 

5 . 4  Method 

There  are  four  major  arrays  that  are  used  to  support 
the  logic  of  this  phase: 

1.  PARAM  (35,10)  =  Project  Underbrush  encoded 
data  handled  in  blocks  of  ten  records  to 
optimize  use  of  memory  end  tape  length. 

2.  NAME  (35)  =  Name  of  each  item  in  the  Project 

Underbrush  encoded  data  as  established  in  a 

Fortran  IV  data  statement. 
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3.  I TAB LF  (5,  35,  10)  ■  Search  parameter  and  inter¬ 
val  card  set.  The  five  represents  one  part 
search  parameter  and  four  parts  interval  in 
order  to  allow  for  the  nine-digit  brightness 
words.  The  ten  allows  for  that  many  case  re¬ 
quests  per  user. 

4.  ITALLY  (3,  10)  ■  Tally  array  where  the  first 
part  holds  a  position  for  each  case  requested 
by  one  user  up  to  ten  positions.  The  second 
part  holds  a  count  of  "and's"  (condition)  for 
each  of  the  above  cases  and  the  third  part 
holds  a  flag  for  each  of  the  cases  to  indicate 
whether  the  interval's  lower  limit  was  found 

to  be  greater  than  its  upper  limit  for  any  "and" 
(condition)  of  a  case,  thus  necessitating  dele¬ 
tion  of  the  entire  case  in  any  subsequent  searches. 

The  logic  proceeds  as  follows: 

Initialize  the  third  (flag)  position  of  ITALLY 
array. 

Read  the  retrieval  request  card  (IPUNCH,  IKEEP) . 

Read  the  search  parameters  and  intervals  into 
ITABLE  array  in  "or"  blocks  (cases)  by  the 
temporary  storage  on  tape  of  each  interval  while 
determining  it  to  be  integer  or  alphanumeric  data. 
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The  item  is  then  read  back  from  tape  and  printed 
in  format. 

Maintain  parts  one  and  two  of  TALLY  array. 

Rewind  scratch  tape. 

Initialize  count  of  records  on  master  Project 
"Underbrush"  (P.U.)  data  tape. 

Read  a  block  (ten  records)  of  P.U.  data  from 
tape . 

Test  to  see  if  it  is  t he  last  block  on  tape. 

Test  the  first  part  of  ITABLE  record  against 
all  of  NAME  array.  When  a  match  is  found,  the 
index  of  the  match  position  allows  the  corres¬ 
ponding  index  position  to  be  located  in  PARAM 
array  record  in  question. 

The  remaining  four  parts  of  ITABLE  record 
(interval)  are  tested  for  content  of  PARAM  item 
of  record  in  question.  If  it  is  satisfied, 
return  to  ITABLE  array  for  the  second  record 
and  so  on  until  the  entire  case  has  been  satis¬ 
fied  at  which  time  the  retrieved  record  is  out¬ 
put.  If  any  condition  is  not  satisfied  by 
the  PARAM  record  in  question,  the  entire  case 
is  skipped  and  execution  branches  to  the  first 
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condition  in  the  next  case. 


Calling  Sequence 

CALL  RETRIV  (ITAPE,  1PABC1,  IPABC2,  MTAPE ,  KEEP) 


Where , 


ITAPE  -  the  tape  unit  on  which  the  data  tape 
finally  appears. 

IPABC1*  the  number  of  completely  filled  blocks 
of  data  on  tape. 


IPABC2-  the  number  of  two-card  data  sets  in  the 
unfilled  block. 


MTAPE  -  the  output  tape  of  retrieved  records  and 
case  numbers. 


KEEP  ■  the  number  of  retrieved  records  on  tape 
(MTAPE) . 


Common  Storage 


None . 


Subroutines  Required 


BRITE 


Accuracy 


Single  precision  floating  point  word. 
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5 . 9  Error  Messages 

Two  error  messages  may  result  from  this  subroutine: 

"NO  MATCH  HAS  BEEN  FOUND  FOR  A6  -  CASE  U  SEARCH 
PARAMETER  12  CHECK  SPELLING” 

occurs  when  the  search  parameter  of  a  search  criterion  card  is  in 
the  improper  field  or  spelled  incorrectly. 

"LOWER  LIMIT  GREATER  THAN  UPPER  LIMIT  IN  CASE  U 
SEARCH  PARAMETER  12 , 

THIS  ENTIRE  CASE  WILL  BE  IGNORED" 
occurs  when  limits  of  the  search  parameter  on  a  search  criterion 
card  are  inverted.  The  lower  limit  should  appear  before  the 
upper  limit  on  the  criterion  cards. 

5.10  Input  Format 

Input  is  provided  to  the  routine  through  cards  and  tape. 


The  card  input  format  is  as  follows: 

CARD  1  -  RETRIEVAL  OPTION  CARD 
COLUMN  ITEM  FORMAT 

1  Punch  option  (IPUNCH)  II 

2  Magnetic  tape  (IKEEP)  II 

CARD  SET  1  -  SEARCH  PARAMETER  AND  SEARCH  CRITERION  SET 

COLUMN  ITEM  FORMAT 

1-6  Search  parameter  (INAM)  A6 

7-9  Last  set  flag  (IEND)  13 

Right  justified 
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12  -  21 


1  Lower  limit  of  interval 
Right  justified 
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2  Lower  limit  of  whole 
part  of  word 

24-33  1  Upper  limit  of  interval  110 

Right  justified 

2  Lower  limit  of  decimal 
part  of  word 

36-45  2  Upper  limit  of  whole  part  110 

of  word 

48-57  2  Upper  limit  of  decimal  110 

part  of  word 

Where  the  2  indicates  the  item  description  when  the  intervals  are 
represented  by  words  of  7  character  length  or  more;  such  as, 
standard  brightness  and  background  brightness  each  of  which  con¬ 
tains  9  characters  (four  whole  number  digits  and  5  decimal  digits). 

The  above  card  represents  one  "and' specif ication  of  a 
possible  35  available  for  retrieval  of  each  "Project  Underbrush" 
record  stored  on  magnetic  tape.  For  example,  a  typical  search 
may  involve  location  of  all  entries  containing  solar  radiation 
values  obtained  under  clear  sky  conditions.  The  series  of  cards 
to  describe  such  a  case  are: 

SOLRAD  1  200 

SKYCOV  1  1 

This  card  group  must  be  terminated  by  a  blank  card.  There  may  be 
as  many  as  ten  such  groups  within  each  set.  The  set  must  be 
terminated  by  a  card  containing  characters  as  follows: 


COLUMNS 


ITEM 


FORMAT 


1-6  Terminal  nines,  999999  A6 

9  Last  set  flag  ( IEND)  13 

Right  justified 

IEND  *  0,  there  are  more  sets  of  search  criteria 
to  examine. 

IEND  =  1,  there  are  no  more  sets  of  search  criteria 
to  examine,  return. 

A  scratch  tape  is  needed  for  temporary  storage  of  the  interval 
fields  while  the  search  parameter  name  is  tested  for  specification 
of  alpha  or  integer  data.  If  the  search  parameter  name  is  one  that 
describes  integer  interval  data,  a  single  format  is  used  to  print 
out  that  name  and  its  interval  data.  If  the  search  parameter  name 
is  one  that  describes  alpha  interval  data,  one  of  four  formats  is 
used  to  print  out  that  name  and  its  interval  data.  They  are  as 
follows : 

Write  out  one  alpha  character  search  criteria 
two  alpha  character  search  criteria 
three  alpha  character  search  criteria 
four  alpha  character  search  criteria. 


5.11  Output  Format 

Subroutine  RETRIV  outputs  two  initial  messages: 

"THIS  IS  THE  RETRIEVAL  PROGRAM" 

"MOUNT  DATA  TAPE  ON  UNIT  1" 
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The  search  parameter  and  search  criterion  table  is  printed  case-by¬ 
case.  At  this  point,  if  any  instance  of  lower  interval  limit 
greater  than  upper  limit  occurred  within  the  table,  the  location 
is  cited. 


A  title,  "RETRIEVED  RECORDS,"  is  printed  as  a  page 
heading  with  six  lines  of  132  Hollerith  characters  giving  the  five 
or  six  character  names  describing  each  item  in  the  "Project  Under- 


brush"  data. 

These  names  are  printed  vertically  atop  the  column 

of  data  they  describe.  Before  each 

retrieved 

record  is  printed 

appears  the 

case  number  to  which  the 

retrieval 

applies.  A  retrieved 

record  appears  on  a  single  printed  line  in  the 

following  format: 

ABBREVIATED 

COLUMNS 

ITEM 

NAME 

FORMAT 

2  -  6 

Mission  Name 

MISNAM 

A5 

7  -  10 

Mission  Number 

MISNUM 

A4 

12 

Sensor  Type 

SENTYP 

11 

14  -  15 

Sensor  Designation 

SENDES 

12 

17  -  22 

Date 

I  DATE 

16 

24  -  26 

Aircraft 

AIRCRF 

A3 

28 

Location 

LOCAT 

A1 

30  -  33 

Sunrise 

SUNRIZ 

14 

35  -  38 

Sunset 

SUNSET 

14 

40  -  43 

Moonrise 

MUNRIZ 

14 

45  -  48 

Moonset 

MUNSET 

14 

50  -  53 

Time  (Local) 

LOCTIM 

14 

55 

Sky  Cover 

SKYCOV 

11 

40 


111   'wm 


1 

1 

57  -  58 

Height  (Clouds) 

CLHITE 

12 

f 

60 

Type  (Clouds) 

CLTYPE 

11 

i 

62  -  63 

Wind  Direction 

WINDIR 

A2 

! 

* 

65  -  66 

Wind  Speed  (Knots) 

WINSPD 

A2 

68  -  69 

Weather  Status 

WESTAT 

A2 

- 

71  -  73 

Wet  Bulb  (F) 

WETBUB 

13 

! 

75  -  77 

Dry  Bulb  (F) 

DRYBUB 

13 

i 

79  -  81 

Relative  Humidity  1 

RELHUM 

13 

! 

83  -  85 

Solar  Radiation 

SOLRAD 

13 

87 

Precipitation 

PRECIP 

Ai 

j 

89 

Ground  Visibility 

GRDVIS 

11 

| 

91  -  92 

Dew  Point  (F) 

DEWPTF 

12 

94  -  96 

Reference  Cube 

REFCUB 

13 

i 

i 

98  -  100 

Background  Results 

BGRRES 

13 

i 

1 

102-  104 

Background  Temp  (C) 

BGRTEM 

13 

\ 

\ 

106 

Moon  Phase 

MUNPHA 

11 

[ 

E 

108-  116 

Standard  Brightness 

STBRIT 

19 

118-  126 

Background  Brightness 

BGBRIT 

19 

I 

! 

128-  130 

Sun/Moon  Shadow 

SUMOSH 

13 

f 

132 

Foliage 

FOLIAG 

11 

f 

i 

The 

use  of  the  "punch"  option  allows  the 

cases  of  anyone 

or  all  users 

to  be  punched  in  the  two- 

card  original 

input  format 

to  subroutine  STORE. 


j 


! 
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The  use  of  the  "keep"  option  allows  the  case  number  plus 
record  satisfying  it  to  be  written  onto  magnetic  tape  in  the 
following  format: 
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WORD 


ITEM 


FORMAT 


1 

Case  number  of  retrieved 
record 

12 

2 

Mission  name 

A5 

3 

Mission  number 

A4 

4 

Sensor  type 

11 

5 

Sensor  designation 

12 

6 

Date 

16 

7 

Aircraft 

A3 

8 

Location 

A1 

9 

Sunrise 

14 

10 

Sunset 

14 

11 

Moonrise 

14 

12 

Moonset 

14 

13 

Local  time 

14 

14 

Sky  cover 

11 

15 

Height  (clouds) 

12 

16 

Type  (clouds) 

11 

17 

Wind  direction 

A2 

18 

Wind  speed 

A2 

19 

Weather  status 

A2 

20 

Wet  bulb 

13 

21 

Dry  bulb 

13 

22 

Relative  humidity 

13 

23 

Solar  radiation 

13 

24 

Precipitation 

11 

25 

Ground  visibility 

12 
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26 

Dew  point 

13 

27 

Reference  cube 

13 

28 

Background  results 

13 

29 

Background  temperature 

11 

30 

Moon  phase 

14 

31 

Standard  brightness 

19 

32 

Background  brightness 

19 

33 

Sun/Moon  shadow 

13 

34 

Foliage 

11 

43 


RETRIV 


MTAPE  =  2 
I  TEST*  1PA6C2 


/  WRITE'.  N 
/'THIS  IS  THE 
RETRIEVAL 
\  PROGRAM' 

ZCZI 

/WRITE:  \ 
/'MOUNT  DATA 
tape  on 
\  UNIT  r  i 


>/ WRITE.:  ' 
/SEARCH  FhfiA- 

(ft\6T«  AMO  S£A«H 
NjRtTMlOM  TABIC' 


1 

Hit *VISL  / 

N  * 

\ 

y-J 

/  WRITE 
'CASE  CN)' 


PAUSE 
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RETRIV 


WRITE. ; 
SEA RCM 


PAR/MAtTEE  < 
INTERVAL 


WRITE: 

'search 
^PARAMETER , 
INTERVAL 


•wnrwf^w  iwiwiWiWi 


RETRIV 


/DO  LOOP 

J-i,  I  pox 


/  DO  LOO P\ 

fw*  i,nuwoks; 


NU/WANft.ITAU^ 


/  DO  LOOP  \ 
(iW*  l#  NliAAAND, 


/DO  LOOP 

w*l,  35* 


N£  CITHER 
VBRICIHTv 


IT/MV-E.  \ 
OjlM.lN) 


C*LL  ftRlTE. 


S  WRITE  > 

/sH0  WATCH.,, 
FOUND..., 
V«e<A  spcuimcv 


STOP 


RETRIV 


6. 


BRITE 


nW'm  irsf  yijrr* 


6. 1  Summary 

This  subroutine  handles  interval  tests  on  integer  words 
which  exceed  six  digits  and  have  been  divided  into  two  words  of 
six  or  less  digits  each. 

6 . 2  Restrictions/Assumptions 
None. 

6 . 3  Options 
None . 

6.4  Method 

Essentially  two  tests  are  performed  in  this  subroutine: 

1 .  Limit  Test 

a.  Test  if  lower  limit  whole  number  part  is 
greater  than  the  upper  limit  whole  number 
part. 

b.  Test  if  lower  decimal  part  is  greater  than 
the  upper  limit  decimal  part. 

2 .  Interval  Test 

a.  Test  if  whole  number  part  of  brightness 
word  falls  between  the  limits  of  the  inter¬ 
val  . 

b.  Test  if  decimal  part  of  brightness  word 
falls  between  the  limits  of  the  interval. 


I 

I 
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Call  BRITE 
Where , 

I  STOW 


ISRCH 


IACC 


I  STOW,  ISRCH,  I ACC) 

The  "Project  Underbrush"  brightness  data 
stored  on  tape;  dimensioned  2  where, 

ISTOW  (1)  =  the  whole  number  part  of 
the  word 

ISTOW  (2)  =  the  decimal  part  of  the  word 

The  interval  limits  from  the  search  para¬ 
meter  and  search  criterion  table;  dimen¬ 
sioned  4  where, 

ISRCH  (1)  =  lower  limit  whole  number 
ISRCH  (2)  =  lower  limit  decimal 
ISRCH  (3)  =  upper  limit  whole  number 
ISRCH  (4)  *  upper  limit  decimal 

A  flag  to  return  acceptance  or  rejection  of 
parameter  interval  fit. 

IACC  =  1,  means  that  both  the  whole  num¬ 
ber  and  decimal  parts  of  a  nine 
character  brightness  value  fall 
between  their  respective  search 
criterion  limits  ***  accept. 

IACC  =  0,  means  that  the  above  was  not 
satisfied  ***  no  accept. 

IACC  =  2 ,  means  that  the  whole  number  and/ 
or  the  decimal  part  of  the  search 
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criterion  interval  is  greater  than 
the  whole  number  and/or  decimal 
part  of  the  lower  limit  of  the 
search  criterion  interval  ***  print 
delete  case. 


6  Common  Storage 
None . 

7  Subroutine  Required 
None . 
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7. 


CORREL 


7 . 1  Summary 

This  subroutine  performs  a  correlation  between  any  two 
sets  of  attributes  for  which  correlation  is  meaningful. 

7 . 2  Res  frictions /As  sumptions 

It  is  assumed  that  the  magnetic  tape  input,  generated 
by  the  IKEEP  retrieval  option,  occurs  through  the  retrieval  phase 

7 . 3  Options 
None . 

7 . 4  Method 

This  subroutine  reads  a  correlation  request  card  con¬ 
taining  the  case  number  of  the  retrieved  records  to  be  used  in 
correlation,  the  abscissa,  and  the  ordinate  of  correlation.  The 
name  array  is  then  searched  for  index  positions  of  abscissa  and 
ordinate  and  when  they  are  found  a  magnetic  tape  of  retrieved 
records  is  read.  These  records  are  tested  for  case  number.  As 
each  pertinent  case  number  is  found,  the  abscissa  and  ordinate 
values  are  placed  in  array  ICOORD  in  memory  and  counted.  Sub¬ 
routine  MISR  is  called  to  calculate  the  means,  variances,  product 
moment  correlation  coefficients,  intercepts  (A)  and  regression 
coefficients  corresponding  to  the  values  of  intercepts  contained 
in  the  output  matrix  above  (B)  .  The  equation  Y=A+BX  is  solved 
for  the  regression  line  coordinates  and  plotted  (PLOT) . 
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In  the  event  that  more  correlations  based  on  the  retrie¬ 
val  case  are  desired,  execution  branches  to  the  beginning  of  the 
subroutine,  rewinds  the  tape,  and  repeats  execution. 

7.5  Calling  Sequence:  CALL  CORREL  (MTAPE,  KEEP) 

Where , 

MTAPE  is  the  input  magnetic  tape  of  retrieved 
records  from  RETRIV 

KEEP  is  the  number  of  records  on  MTAPE. 

7 . 6  Common  Storage 
None . 

7 . 7  Subroutines  Required 
ARRAY,  MI SR,  PLOT. 

7 . 8  Accuracy 

Single  precision  floating  point  word. 

7 . 9  Error  Messages 

"YOU  HAVE  PROBABLY  PICKED  A  CASE  NUMBER  WHICH  HAS 
NO  RETRIEVED  RECORD  ON  TAPE" 

is  printed  when  the  entire  tape  has  been  read  and  no  case  request 
has  been  satisfied. 

7.10  Input  Format 

CORRELATION  REQUEST  CARD 

COLUMN  ITEM  FORMAT 

1  Case  number  of  retrieved  rec-  II 

ord  to  be  used  in  correlation 
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7  - 

12  Abscissa 

A6 

14  - 

19  Ordinate 

A6 

25 

0,  there  are  more  cases  to  be  A1 
correlated;  continue  in  the 
correlation  cycle. 

1,  there  are  no  additional 
to  be  correlated;  end  the 
tation. 

cases 

compu- 

RETRIEVED  RECORD  TAPE 

WORD 

ITEM 

FORMAT 

1 

Case  number  of  retrieved 
record 

12 

2 

Mission  name 

A5 

3 

Mission  number 

A4 

4 

Sensor  type 

11 

5 

Sensor  designation 

12 

6 

Date 

16 

7 

Aircraft 

A3 

8 

Location 

A1 

9 

Sunrise 

14 

10 

Sunset 

14 

11 

Moonrise 

14 

12 

Moons et 

14 

13 

Local  time 

14 

14 

Sky  cover 

11 

15 

Height  (clouds) 

12 

16 

Type  (clouds) 

11 

17 

Wind  direction 

A2 
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18 

Wind  speed 

A2 

19 

Weather  status 

A2 

20 

Wet  bulb 

13 

21 

Dry  bulb 

13 

22 

Relative  humidity 

13 

23 

Solar  radiation 

13 

24 

Precipitation 

11 

25 

Ground  visibility 

12 

26 

Dew  point 

13 

27 

Reference  cube 

13 

28 

Background  results 

13 

29 

Background  temperature 

11 

30 

Moon  phase 

14 

31 

Standard  brightness 

19 

32 

Background  brightness 

19 

33 

Sun/moon  shadow 

13 

34 

Foliage 

11 

Output  Format 

This  subroutine  outputs  a  computer  listing  containing: 

1.  Correlation  request  card  information 

2.  Number  of  retrieved  records  on  tape 

3.  X-mean  and  Y-mean 

4.  X- variance  and  Y-variance 

5.  Product -moment  correlation  coefficients 

6.  Table  of  regression  line  coordinates 

7.  Plot  of  the  regression  line. 


I 
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8 .  ARRAY 


8 . 1  Summary 

This  subroutine  converts  data  from  single  to  double 
dimension  or  vice  versa.  The  program  is  used  to  link  the  user 
program  which  contains  double  dimension  arrays  and  other  subrou¬ 
tines  which  operate  on  arrays  of  data  in  a  vector  fashion.  This 
program  was  provided  by  the  IBM  program  library. 


8. 2  Restrict ions/ Assumptions 

Inputs  to  this  subroutine  occur  through  magnetic  tape 
from  the  retrieval  phase. 

8 . 3  Options 
None . 


8.4  Calling  Sequence:  CALL  ARRAY  CMODE ,  I,  J,  N,  Mt  S,  D) 

Where , 

MODE  =  Code  indicating  type  of  conversion. 

1  =  From  single  to  double  dimension. 

2  =  From  double  to  single  dimension. 

I  =  Number  of  rows  in  actual  data  matrix. 

J  =  Number  of  columns  in  actual  data  matrix. 

N  =  Number  of  rows  specified  for  the  matrix  D 
in  the  dimension  statement. 

M  =  Number  of  columns  specified  for  the  matrix  D 
in  the  dimension  statement. 
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S  ■  If  mode  =  1,  this  vector  is  input  which  con¬ 
tains  the  elements  of  a  data  matrix  of  size 
I  by  J. 

Column  I  +  1  of  data  matrix  follows  Column 
I,  etc. 

If  mode  =  2,  this  vector  is  output  repre¬ 
senting  a  data  matrix  of  size  I  by  J  con¬ 
taining  its  columns  consecutively. 

D  =  If  mode  =  1,  this  matrix  of  size  N  by  M  is 
output,  containing  a  data  matrix  of  size  I 
by  J  in  the  first  I  rows  and  J  columns. 

If  mode  =  2,  this  N  by  M  matrix  is  input  con¬ 
taining  a  data  matrix  of  size  I  by  J  in  the 
first  I  rows  and  J  columns. 

Common  Storage 
None. 

Subroutines  Required 
None . 

Accuracy 

Single  precision  floating  point  word. 

Error  Messages 
None. 


64 


9. 


MISR 


9.1  Summary 

This  subroutine  computes  means,  standard  deviation, 
skewness  and  kurtosis,  correlation  coefficients,  regression  coef¬ 
ficients,  and  standard  errors  of  regression  coefficients  when 
there  are  missing  data  prints.  The  user  identifies  missing  data 
by  means  of  a  numeric  code.  Those  values  having  this  code  are 
skipped  in  computing  the  statistics.  In  the  case  of  the  correla¬ 
tion  coefficients,  any  pair  of  values  is  skipped  if  either  one  of 
them  ai :  missing.  This  program  was  provided  by  the  IBM  program 
library. 

9. 2  Restrict  ions/ As sumptions 

The  retrieval  phase  must  be  executed  before  this  sub¬ 
routine  can  be  used.  Also,  this  subroutine  cannot  distinguish  a 
blank  and  a  zero.  Therefore,  if  a  blank  is  specified  as  a  missing 
data  code  in  input  cards,  it  will  be  treated  as  a  zero. 

9.3  Options 

None . 


9.4  Method 

Least  squares  regression  lines  and  product -moment  cor¬ 
relation  coefficients  are  computed. 


9.5  Calling  Sequence 

Call  MISR  (NO,  M,  X,  CODE,  XBAR,  STD,  VAR,  SKEW,  CURT, 
R,  N,  A,  B,  S,  ILR)  where, 
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NO  ■  Number  of  observations. 

M  *  Number  of  variables. 

X  -  Input  data  matrix  of  size  NO  by  M. 

CODE  ■  Input  vector  of  length  M,  which  contains  a 
numeric  missing  data  code  for  each  variable. 
Any  observation  for  a  given  variable  having 
a  value  equal  to  the  code  will  be  dropped. 

XBAR  *  Output  vector  of  length  M  containing  means. 

STD  =*  Output  vector  of  length  M  containing  standard 

deviations . 

VAR  =*  Output  vector  of  length  M  containing  vari¬ 
ance  . 

SKEW  =  Output  vector  of  length  M  containing  skew¬ 
ness. 

CURT  =  Output  vector  of  length  M  containing  kurto- 
sis . 

R  =  Output  matrix  of  product -moment  correlation 
coefficients.  This  will  be  the  upper  tri¬ 
angular  matrix  only,  since  the  M  by  M  matrix 
of  coefficients  is  symmetric.  (Storage  Mode  1) 
N  =  Output  matrix  of  number  of  pairs  of  observa¬ 
tions  used  in  computing  the  correlation  coef¬ 
ficients.  Only  the  upper  triangular  portion 
of  the  matrix  is  given.  (Storage  Mode  1) 

A  =  Output  matrix  (M  by  M)  containing  intercepts 
of  the  regression  lines  (A)  of  the  form 

Y=A+BX.  The  first  subscript  of  this  matrix 
67 


refers  to  the  independent  variable  and  the 
second  to  the  dependent  variable.  Matrix  A 
is  stored  in  a  vector  form. 

B  *  Output  matrix  (M  by  M)  containing  regression 
coefficients  (B)  corresponding  to  the  values 
of  intercepts  contained  in  the  output  ma¬ 
trix  A. 

S  =  Output  matrix  (M  by  M)  containing  standard 
errors  of  regression  coefficients  corres¬ 
ponding  to  the  coefficients  contained  in 
the  output  matrix  B. 

IER  =  0,  no  error. 

1,  if  number  of  non-missing  data  elements 
for  the  Jth  variable  is  two  or  less.  In 
this  case,  STD(J) ,  SKEW(J) ,  and  CURT(J)  are 
set  to  10  **  75.  All  values  of  R,  A,  B, 
and  S  related  to  this  variable  are  also  set 
to  10  **  75. 

2,  if  variance  of  Jth  variable  is  less  than 
10  **  (-20).  In  this  case,  STD(J),  SKEW(J) , 
and  CURT(J)  are  set  to  10  **  75.  All  values 
of  R,  A,  B,  and  S  related  to  this  variable 
are  also  set  to  10  **  75. 


Common  Storage 
None . 
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Subroutines  Required 
None . 

Accuracy 

Single  precision  floating  point  word. 

Error  Message 
None . 


69 


• *  *  Mtur/iHaii— 


10. 


PLOT 


10.1 


Summary 


This  subroutine,  provided  by  the  IBM  program  library, 


is  used  for  plotting  on  a  printer.  It  will  plot  a  graph  with  one 
independent  variable  and  up  to  nine  dependent  variables,  with  the 
additional  ability  to  plot  a  calculated  curve. 


10.2 


Res t rict ions /As sump t ions 

The  retrieval  phase  must  be  executed  before  this  sub¬ 


routine  can  be  used. 


10 . 3  Options 

None . 

10 . 4  Method 

The  program  first  prints  the  chart  number.  It  then  cal¬ 
culates  the  extremes  of  X  and  Y,  unless  these  are  specified.  A 
rough  scale  is  calculated,  assuming  50  divisions  vertically  and 
100  horizontally. 

The  print  position,  YPR,  is  set  to  the  maximum  value  of 
Y,  and  all  independent  variable  sets,  including  calculated  values, 
are  searched  for  values  of  Y  within  one-half  scale  division  of  YPR. 
If  any  are  found,  an  integer  (1-9)  is  printed  on  that  line  at  the 
proper  X  value.  The  value  of  the  integer  is  that  of  the  set  of 
dependent  variables  hose  Y  value  was  matched.  \n  asterisk  is 
printed  for  a  calculated  value.  The  search  is  done  first  on  cal¬ 
culated  value,  then  on  ascending  sets  of  Y.  Only  the  last  match 


70 


found,  if  any,  prints. 


The  print  position  of  YPP.  is  decremented  by  the  scale, 
and  the  process  repeated  until  the  minimum  value  of  Y  is  reached. 

10.5  Calling  Sequence 

Call  PLOT  (NO,  A,  N,  M,  NFUNC,  FUNC,  XLAX,  XLIN,  YLAX, 
YLIN)  where, 

NO  =  A  fixed  point  number,  up  to  three  digits, 
printed  as  the  chart  number. 

A  *  A  vector  whose  first  N  positions  contain 
the  independent  variable,  and  whose  next 
M  sets  of  N  positions  contain  the  dependent 
variables . 

N  =  The  number  of  observations  for  each  of  the 
variables.  If  N“0,  this  means  no  tabulated 
points  are  to  be  plotted. 

M  =  The  number  of  variables  (dependent  and 
independent) . 

NFUNC  =  Zero  or  negative  if  no  calculated  curve 
wanted,  positive  if  one  is  wanted. 

FUNC  =  Subroutine  of  the  form  SUBROUTINE  FUNC  (X,Y) 
which  returns  a  value  of  Y  when  given  a 
value  of  X.  The  program  calling  PLOT  must 
have  an  EXTERNAL  FUNC  statement. 

XLAX,  XLIN,  YLAX,  YLIN  ■  Maximum  and  minimum  values 
of  the  independent  and  dependent  variables 


to  be  used  in  the  plot  if  XLAX-XLIN.  The 
program  calculates  its  own  maximum  and 
minimum  for  the  independent  variable. 
Program  operates  similarly  for  YLAX=YLIN. 

Common  Storage 
None. 


Subroutines  Required 
FUNC  (if  used)  and  SCAL. 

Accuracy 

Single  precision  floating  point  word. 
Error  Message 
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11  •  SCAI, 

11.1  Summary 

This  final  subroutine,  an  integral  part  of  the  PLOT  sub¬ 
routine,  is  used  to  reduce  t'ne  PLOT  scale  l»y  roundjng  the  values. 

11.2  tfe  strictions /hj_  sumptions 

This  program  ir-  used  in  conjunction  vith  the  PICT  routine. 


\ 

11.3  Options 
None. 


11.4  Method 

The  rough  scale  calculated  b?  t)  >  TLO!  '•culi.  e  is  reduced 
to  l,  2,  2.5,  or  5  times  10  to  an  integral  pown .  The  end  points 
are  adjusted  to  be  multiples  of  'he  seme. 


'•  v' 


12. 


DATA  REDUCTION 


Part  of  the  UNDERBRUSH  Test  Range  support  activity  con¬ 
cerns  the  collecting  and  recording  of  ground  truth  information  in 
conjunction  with  sensor  flight  test  operations.  The  type  and  ex¬ 
tent  of  the  data  depend  primarily  upon  the  nature  of  the  sensor 
system  being  tested  and  the  level  of  detail  needed  to  support 
accurate  evaluation.  The  ground  truth  information  obtained  can 
be  divided  into  three  segments:  mission  data,  target  information, 
and  atmospheric,  observations .  Mission  data  pertain  essentially 
to  aircraft  flight  information  such  as  heading,  altitude,  time 
over  target,  and  the  start  and  end  times  for  each  sortie.  Target 
information  typically  includes  ihe  type,  location,  orientation, 
and  dimensions  of  each  target  and,  when  required,  such  specialized 
information  as  surface  reflectivity  and  thermal  measurements.  Data 
on  cloud  condition,  temperature,  precipitation,  wind,  and  solar 
radiation  are  included  within  the  atmospheric  observation  category. 

The  objective  of  the  Computer  Assisted  Target  Analysis 
program  was  to  design  and  test  a  computer  program  to  store,  retrieve, 
and  correlate  selected  ground  truth  information.  Before  the  compu¬ 
ter  program  could  be  created,  however,  selection  of  ground  data 
parameters  had  to  be  accomplished.  Unfortunately  the  types  of 
"round  information  being  collected  and  the  manner  in  which  they  are 
presented  has  not  remained  constant.  Hence,  one  of  the  first  tasks 
in  the  program  was  to  search  through  the  UNDERBRUSH  ground  truth 
files  and  list  all  the  information  categories.  From  this  list  32 
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data  parameters  were  selected  for  inclusion.  A  short  description 
of  these  parameters  follows. 

1.  Mission  Name  and  Number  -  This  nomenclature  is 

included  only  if  the  imagery  is  on  file  at 
the  RADC/EMIRC  Reconnaissance  Data  Base. 

2.  Sensor  Type  -  Type  of  sensor,  e.g.,  Radar,  IR, 

Photo,  Laser. 

3.  Sensor  Designation  -  Sensor  system  identifica¬ 

tion  name. 

4.  Date  -  Date  of  ground  truth  observation. 

5.  Aircraft  -  Type  of  aircraft  carrying  sensor 

equipment. 

6.  Location  -  Site  where  atmospheric  observations 

were  collected.  (See  Appendix  A) 

7.  Sunrise  -  Local  time  of  sunrise. 

8.  Sunset  -  Local  time  of  sunset. 

9.  Moonrise  -  Local  time  of  moonrise. 

10.  Moonset  -  Local  time  of  moonset. 

11.  Time  -  Local  time  of  each  recorded  observation. 

12.  Sky  Cover  -  Amount  of  cloud  cover. (See  Appendix  A) 

13.  Cloud  He\g  t  -  Estimated  height  of  clouds  above 

the  ground  in  thousands  of  feet. 

14.  Cloud  Typ'3  -  Type  of  cloud  formation. 

(See  Appendix  A) 

15.  Wind  Direction  -  Wind  direction  relative  to  the 

eight  points  of  the  compass. 


16.  Wind  Speed  -  Velocity  in  knots. 

17.  Weather  Status  *  Weather  phenomena  at  time  of 

observation.  (See  Appendix  A) 

18.  Wet  Bulb  -  Wet  bulb  temperature  reading  in 

tenths  of  degrees  F. 

19.  Dry  Bulb  -  Dry  bulb  temperature  reading  in 

tenths  of  degrees  F. 

20.  Relative  Humidity  -  Percent  of  relative  humidity. 

21.  Solar  Radiation  -  Reading  of  solar  radiation 

recorded  in  hundredths. 

22.  Precipitation  -  Amount  of  precipitation  occurring 

between  the  time  of  observations. 

23.  Ground  Visibility  -  Estimate  of  visibility. 

24.  Dew  Point  -  Dew  point  temperature  in  tenths  of 

degrees  F. 

25.  Reference  Cube  -  Reference  temperature  value  in 

tenths  of  degrees  C  applied  to  radiometric 
readings . 

26.  Background  Results  -  Background  reading  of 

radiometer  in  tenths  of  degrees  C. 

27.  Background  Temperature  -  Background  temperature 

in  tenths  of  degrees  C. 

28.  Moon  Phase  -  Moon  phase  to  nearest  quarter. 

29.  Standard  Brightness  -  Reference  value  applied 

to  footlamberts  reading. 
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30.  Background  Brightness  •  Background  brightness 

in  footlamberts . 

31.  Sun/Moon  Shadow  -  Direction  of  shadow  angle  in 

degrees . 

32.  Foliage  -  Amount  of  foliage  coverage  at  place 

of  recorded  observation.  (See  Appendix  A) 

In  order  to  minimize  the  amount  of  computer  storage,  yet 
include  complete  data,  alphabetic  and  numeric  codes  were  developed. 
Appendix  A  lists  these  codes  along  with  their  representative  valuer 
The  encoded  ground  information  was  transcribed  onto  a  format  sheet 
and  then  to  a  keypunch  data  card.  Since  two  cards  were  needed  to 
record  a  set  of  data  for  a  given  observation  time,  a  sequential 
number  was  assigned  to  each  punch  card  to  ensure  correct  order. 

The  punch  card  format  is  as  follows: 

CARD  1 


COLUMNS 

ITEM 

FORMAT 

1 

-  5 

MISNAM 

-  Mission  Name 

AS 

6 

-  9 

MISNUM 

■  Mission  Number 

A4 

11 

SENTYP 

-  Sensor  Type 

11 

13 

-  14 

SENDES 

-  Sensor  Designation 

12 

16 

-  21 

I  DATE 

-  Date 

16 

23 

-  25 

AIRCRF 

-  Aircraft 

A3 

27 

LOCAT 

-  Location 

A1 

29 

-  32 

SUNRIZ 

-  Sunrise 

14 

34 

-  37 

SUNSET 

•  Sunset 

14 

77 


39 

-  42 

MUNRIZ 

-  Moonrist* 

14 

44 

-  47 

MUNSF.T 

-  Moonset 

14 

49 

-  52 

LOCTIM 

-  Time  (Local) 

14 

54 

SKYCOV 

-  Skycover 

1 1 

56 

-  57 

CLHITE 

-  Height  (Clouds) 

12 

59 

CLTYPE 

-  Type  (Clouds) 

11 

61 

•  62 

WINDIR 

-  Wind  Direction 

A2 

64 

-  65 

WINSPD 

-  Wind  Speed  (Knots) 

A2 

67 

-  68 

WESTAT 

-  Weather  Status 

A2 

70 

-  72 

WETBUB 

-  Wet  Bulb  (F) 

13 

••  77 

-  80 

Sequential  Card  Identification-  14 
Never  Read,  Written,  or  Processed 

CARD  2 

COLUMNS  ITEM  FORMAT 


1 

-  3 

DRYBUB 

-  Dry  Bulb  (F) 

13 

5 

-  7 

RELHUM 

-  Relative  Humidity  \ 

13 

9 

-  11 

SOLRAD 

-  Solar  Radiation 

13 

13 

PRECIP 

-  Precipitation 

A1 

15 

GRDV1S 

-  Ground  Visibility 

11 

17 

-  18 

DEWPTF 

-  Dew  Point  (F) 

12 

20 

-  22 

REFCUB 

-  Reference  Cube 

13 

24 

-  26 

BGRRES 

-  Background  Results 

13 

28 

-  30 

BGRTEM 

-  Background  Temp  (C) 

13 

32 

MUNPHA 

-  Moon  Phase 

1 1 

34 

-  37 

SBRITW 

•  Standard  Brightness 
(Whole) 

14 

78 


38  -  42 


SBRITD  -  Standard  Brightness 
(Decimal) 


15 


44 

-  47 

BBRITW 

-  Background  Brightness 
(Whole) 

14 

48 

-  52 

BBRITD 

-  Background  Brightness 
(Decimal) 

IS 

54 

-  56 

SUMOSH 

-  Sun/Moon  Shadow 

13 

58 

FOLIAG 

-  Foliage 

11 

*«  77 

-  80 

Sequential  Card  Identification-  14 
Never  Read,  Written,  or  Processed 

A  total  of  3,108  observations  were  encoded  and  6,216 
data  cards  were  keypunched,  machine  verified,  sorted  into  proper 
sequence,  and  stored  on  magnetic  tape.  In  addition,  ground  truth 
information  covering  Fort  Sill  and  Fort  Knox  were  also  processed, 
although  those  data  could  not  be  incorporated  as  part  of  the  UNDER 
BRUSH  test  records. 


13.  COMPUTER  TEST  RESULTS 

This  section  presents  the  results  of  the  various  tests 
undertaken  in  conjunction  with  the  RADC  Program  Monitor  to  check 
out  the  operation  of  the  computer  program  and  to  demonstrate  the 
utility  of  the  program  in  support  of  UNDERBRUSH  evaluations.  All 
the  computer  operations  were  conducted  using  the  GE-645  computer 
at  Rome  Air  Development  Center.  For  complete  instructions  con¬ 
cerning  the  set-up  of  the  card  deck  and  the  operational  procedures 
refer  to  Section  2  of  this  report. 

Although  the  following  examples  cover  different  situ¬ 
ations,  the  basic  computer  procedure  remains  the  same.  First  the 
requestor  specifies  the  categories  and  limits  of  information  to  be 
retrieved.  The  computer  *.hen  searches  the  master  magnetic  tape 
containing  all  the  UNDERBRUSH  records  and  retrieves  those  sets  of 
data  which  comply  with  the  requested  parameters.  Retrieval  is 
accomplished  by  copying  the  desired  sets  of  observations  onto  a 
spare  magnetic  tape.  This  spare  tape  is  then  used  to  generate  a 
listing  of  requested  information  as  well  as  the  input  to  the  cor¬ 
relation  program. 

It  should  be  noted  that  because  of  the  nature  of  the 
data  being  collected  for  any  given  set  of  UNDERBRUSH  observatiors 
some  of  the  categories  contain  no  information;  hence,  a  blank  is 
recorded  on  the  master  tape.  However,  the  GE-645  computer  treats 
a  blank  as  a  zero.  Therefore,  during  the  correlation  phase  if 
cither  of  the  parameters  being  analyzed  contains  blanks  the 
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resultant  computation  is  distorted.  To  prevent  this  possibility, 
the  retrieval  case  should  include  the  correlation  perimeters. 

This  ensures  that  the  input  to  the  correlation  phase  contains  com¬ 
plete  data. 

13.1  Dry  Bulb  and  Wet  Bulb  Temperature  Vs.  Solar  Radiation 

The  purpose  of  this  test  was  to  determine  whether  solar 
radiation  is  related  to  temperature,  and  if  so,  whether  the  dry 
bulb  temperature  or  wet  bulb  temperature  is  more  important.  The 
first  step  was  to  retrieve  all  the  UNDERBRUSH  observations  con¬ 
taining  solar  radiation  readings,  then  compare  the  temperature 
readings  with  values  for  solar  radiation.  A  total  of  920  observa¬ 
tions,  S52  dry  bulb  and  368  wet  bulb,  were  retrieved  which  when 
correlated  produced  the  following  results: 


Average  dry  bulb  temperature  70.1*F 

Average  wet  bulb  temperature  59.8#F 

Average  solar  radiation  from 
dry  bulb  observations  0.59 

Average  solar  radiation  from 
wet  bulb  observations  0.63 

Correlation  coefficient  with 
dry  bulb  0. 28 

Correlation  coefficient  with 
wet  bulb  0.03 


Analysis  shows  that  there  is  a  slight  tendency  for  solar  radiation 
to  be  related  to  dry  bulb  temperature  at  least  in  comparison  with 
the  wet  bulb  readings.  However,  it  is  also  evident  from  the  low 
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correlation  coefficient  values  that  ambient  air  temperature  does 
not  vary  proportionately  with  the  intensity  of  solar  radiation. 

13.2  Dry  Bulb  Temperature  Vs.  Solar  Radiation  for  Various 

Sky  Conditions 

In  the  previous  example  solar  radiation  and  temperature 
were  correlated  with  mixed  results.  Here  the  problem  has  been 
further  refined  by  comparing  data  obtained  under  similar  sky 
conditions.  Four  types  of  sky  cover  were  differentiated:  clear, 
scattered,  broken,  and  overcast  conditions. 

Results  of  the  computer  calculations  show  a  distinct 
relationship  between  the  amount  of  cloud  cover  and  solar  radiation 
intensity.  As  would  be  expected  radiation  level  decreases  as  the 
cloud  cover  increases.  By  the  same  token  the  product-moment 
correlation  coefficient  increases  steadily  indicating  that  under 
cloudy  conditions  a  relationship  seems  to  exist  between  radiation 
intensity  and  dry  bulb  temperature.  This  relationship,  however, 
is  not  strong  enough  to  be  causal,  but  rather  may  be  associated 


with  other  atmospheric  factors. 

SKY 

DRY 

SOLAR 

CORRELATION 

CONDITION 

BULB  °F 

RADIATION 

COEFFICIENT 

Clear 

68.1 

0.69 

0.16 

Scattered 

72.4 

0.68 

0.28 

Broken 

76.7 

0.53 

0.42 

Overcast 

64.6 

0.27 

0.41 
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13 . 3  Background  Temperature  Vs.  Background  Radiometer  Results 

The  intensity  of  infrared  radiation  from  a  given  object 

is  a  function  of  its  surface  temperature  and  emissivity.  Since 
the  background  radiation  readings  collected  at  the  UNDERBRUSH  Test 
Range  were  obtained  by  sampling  the  ground,  it  can  be  assumed  that 
the  emissivity  effect  would  be  relatively  constant  and  hence  would 
not  be  much  of  a  factor  in  determining  changes  in  radiation  levels. 
If  the  above  assumptions  are  correct  a  close  correlation  should 
exist  between  background  temperature  and  background  radiation.  To 
chock  the  validity,  all  the  UNDERBRUSH  temperature  and  radiation 
data  were  retrieved  on  tape  and  then  processed  through  the  cor¬ 
relation  routine  to  obtain  the  correlation  coefficient.  The  compu¬ 
ter  results  showed  that  the  product -moment  correlation  coefficient 
was  0.91  indicating  a  very  high  correlation  between  the  background 
temperature  readings  and  background  radiation  levels. 

13.4  Windspeed  Vs.  Background  Radiometric  Readings 

The  objective  of  this  test  was  to  investigate  whether 
the  velocity  of  wind  affects  the  intensity  of  terrestrial  radiation. 
The  first  step  was  to  determine  what  type  of  data  should  be  re¬ 
trieved  and  define  the  parametric  limits.  Since  terrestrial  radi¬ 
ation  is  primarily  a  function  of  surface  temperature,  which  in 
turn  is  strongly  influenced  by  sunlight,  only  those  observations 
recorded  at  night  between  2030  and  2400  hours  were  considered. 

In  addition,  the  data  being  retrieved  were  separated  into  groups 
based  upon  dry  bulb  temperature  readings. 
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Statistical  results  of  the  173  observations  retrieved 


are  as  follows: 


Mid-Point  Dry  Bulb  Temperature 

32 . 5  °F 

47 . 5°F 

62 . 5°F 

77. 5°F 

Average  Radiometric  Reading 

4 . 9°C 

10 . 2°C 

19 . 6°C 

26. 9°C 

Correlation  Coefficient 

0.37 

O 

o 

1 

0.21 

1 

o 

• 

Judging  from  the  rather  erratic  product-moment  correlation  coef¬ 
ficient  values  it  would  appear  that  the  effect  of  wind  on  the 
intensity  of  terrestrial  radiation  is  minor  compared  to  other  fac¬ 
tors.  Indeed  the  increase  of  radiometric  readings  with  respect 
to  increased  air  temperature  is  vividly  demonstrated  by  the  re¬ 
sults.  However,  two-thirds  of  the  radiometric  readings  were 
collected  under  calm  wind  conditions  and  in  only  13  instances  was 
the  wind  speed  greater  than  5  knots.  More  data  are  needed, 
especially  under  higher  wind  velocities  before  conclusive  results 
can  be  obtained. 

13.5  Variable  Relative  Humidity  Vs.  Solar  Radiation 

The  purpose  of  this  final  example  was  to  determine  if 
the  UNDERBRUSH  data  showed  any  correlation  between  the  intensity 
of  solar  radiation  and  relative  humidity.  Since  the  intensity  of 
solar  radiation  varies  throughout  the  daytime,  generally  being 
highest  around  mid-day  and  lowest  at  dawn  and  dusk,  only  those 
solar  radiation  readings  obtained  between  1000  and  1500  hours 
local  time  were  considered.  The  search  parameters  were  further 
refined  by  requesting  that  the  retrieved  data  be  separated  into 
two  groups,  one  having  relative  humidity  values  between  40  and  60 


percent,  and  the  other  between  90  ai»l  i00  percent.  Af the 
search  function  was  completed  83  records  *ere  reciivYed,  70  con¬ 
taining  low  relative  humidity  values  and  only  vitb  high  humidity. 

Results  of  the  correlation  showed  that  for  t  first 
group  the  average  relative  humidity  was  49  percent  and  the  solar 
radiation  value  0.81.  The  second  group  had  an  average  relative 
humidity  reading  of  98  percent  and  a  corresponding  solar  radiation 
value  of  only  0.51.  Thus,  a  solar  radiation  difference  of  0.30 
exists  between  the  two  groups.  The  product -moment  correlation 
coefficient  value  for  the  first  group  was  computed  to  be  -0.19 
indicating  a  slight  inverse  correlation.  Yet  the  second  group 
showed  a  very  strong  inverse  correlation  of  -0.70.  From  the 
data  it  appears  that  although  low  humidity  does  not  influence 
the  intensity  of  solar  radiation  very  much,  high  humidity  is  in¬ 
deed  an  important  factor.  As  the  humidity  level  approaches  the 
saturation  point  the  amount  of  incoming  solar  radiation  decreases. 
This  indicates  that  atmospheric  absorption  of  short  wave  solar 
radiation  is  more  pronounced  as  the  percent  of  water  vapor  in  the 
atmosphere  reaches  a  maximum. 

Unfortunately  the  limited  amount  of  data  available  con¬ 
cerning  high  humidity  values  prohibits  conclusive  analysis,  especi¬ 
ally  since  an  apparently  abnormal  radiatio.i  reading  was  included. 
For  instance,  portions  of  the  five  observations  are  as  follows: 


DATE 

TIME 

RELATIVE 
HUMIDITY C%) 

SOLAR 

RADIATION 

660927 

1330 

100 

0.20 

660927 

1400 

100 

0.25 

660927 

1430 

100 

0.32 

670708 

1330 

96 

0.39 

670708 

1400 

96 

1.40 

As  can  be  seen,  four  of  the  five  solar  radiation  values  are  reason¬ 
ably  close  together,  yet  the  fifth  one  is  way  out  of  line,  especi¬ 
ally  considering  that  the  last  two  observations  were  sequentially 
recorded  one-half  hour  apart  on  the  same  day.  If  the  last  obser¬ 
vation  were  deleted,  the  inverse  correlation  between  solar  radiation 
and  relative  humidity  would  become  even  stronger. 


14. 


CONCLUSIONS  AND  RECOMMENDATIONS 


The  objective  of  the  Computer  Assisted  Target  Analysis 
program,  that  is,  to  design  and  test  a  computer  program  to  store, 
retrieve,  and  correlate  selected  ground  truth  information,  was 
accomplished.  The  UNDERBRUSH  data  gathered  from  the  advent  of  the 
test  range  until  April  1969  were  encoded,  keypunched,  and  stored 
on  magnetic  tape.  In  addition,  atmospheric  observations  covering 
Fort  Sill  and  Fort  Knox  also  were  processed. 

In  the  preceeding  section  specific  examples  tost: 
the  search,  retrieval,  and  correlation  techniques  of  : he  compute 
program  were  described.  The  program  utility  is  by  no  means  limited 
to  those  types  of  data  manipulation.  Some  other  considerations  are 
described  below. 

1.  Since  the  master  tape  contains  data  covering 
three  years  of  UNDERBRUSH  observations,  some 
climatic  information  can  be  obtained  from  the 
program.  For  instance,  one  of  the  problems 
encountered  in  planning  a  sensor  test  exercise 
is  to  estimate  the  number  of  flying  days.  Al¬ 
though  the  UNDERBRUSH  data  were  collected  only 
during  active  operations,  the  amount  of  cloud 
cover  and  ceiling  height  for  a  given  period  can 
be  determined,  thus  providing  an  indication  of 
the  typical  sky  condition  which  can  be  antici¬ 
pated. 


Although  there  are  other  program  functions, 
such  as  storing  additional  data  on  the  master 
tape,  the  search  and  retrieval  function  is 
undoubtedly  the  most  important.  Through  the 
use  of  search  input  cards  the  program  user  can 
define  the  data  to  be  searched  and  retrieved. 
The  actual  retrieval  can  be  output  in  several 
ways:  printed  out  by  the  computer,  transferred 
to  another  magnetic  tape,  or  punched  on  data 
cards.  These  output  options  provide  a  versa* 
tility  which  can  be  useful  for  other  projects. 
An  example  might  be  a  study  on  low  light  level 
TV  operating  requirements  which  could  incor¬ 
porate  the  UNDERBRUSH  Test  Range  nocturnal 
illumination  readings.  The  illumination  data 
could  be  punched  out  on  cards  by  the  computer, 
then  compiled  with  data  from  other  sources. 

The  basic  computer  program  can  be  exploited  for 
other  uses  besides  UNDERBRUSH  d&ta  analysis. 

By  modifying  the  data  categories  a  search  and 
retrieval  system  could  be  utilized  for  such 
projects  as  the  data  base  imagery  file,  multi¬ 
sensor  target  inventory,  and  SEA  support  data 
analysis.  Any  information  that  lends  itself  to 
encoding  can  be  readily  handled. 
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Data  reduction  for  the  Computer  Assisted  Target  Analysis 
progrum  consisted  primarily  of  atmospheric  observations;  however, 
consideration  should  be  given  toward  expanding  the  project  to  in¬ 
clude  specific  flight  and  target  information.  Data  have  already 
been  collected  at  the  UNDERBRUSH  Test  Range  covering  such  diverse 
items  as  aircraft  heading  and  altitude,  sensor  mode  stations, 
radiometric  and  thermal  readings  of  specific  targets,  MTI  informa¬ 
tion,  and  location  and  orientation  of  static  arrays.  The  useful¬ 
ness  of  such  data  for  sensor  evaluations  would  be  enhanced  if  they 
were  more  readily  accessible. 


APPENDIX  A 


A. 

COMPUTER 

ASSISTED  TARGET  ANALYSIS 

A.  1 

Sensor  Type  -  SENTYP 

1  - 

Radar 

2  - 

Infrared 

3  - 

Photography 

4  - 

Laser 

A. 2 

Sensor  Designation  -  SENDES 

Refer  to 

RADC/EMIRC  Reconnaissance 

codes. 

A. 3 

Locat ion 

-  LOCAT 

A  - 

Site  1 

K 

B  - 

Site  2 

L 

C  - 

Site  3 

M 

D  - 

Site  4 

N 

E  - 

SAM  Site 

0 

F  - 

Sampan 

P 

G  - 

General  UNDERBRUSH 

Area 

Q 

H  - 

Triangle 

R 

I  - 

Intercoastal  Waterway 

S 

J  - 

Route  287 

T 

A. 4 

Sky  Cover  -  SKYCOV 

1  - 

Clear 

4 

2  - 

Scattered  Clouds 

5 

3  - 

Broken  Clouds 

ENCODING  KEY 


Data  Base  sensor 

-  POL 

-  Range  70 

-  Holley  Field 

-  Range  C-52 

•  Rock  Hill  Assault 
Strip 

-  Clay  Pit 

-  Range  C-74 

-  Fort  Knox 

-  Range  B-75 

-  Fort  Sill 

-  Overcast 

-  Obscured 
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A.  S 


A. 6 


A.  7 


A.  8 


mi 


Cloud  Type  ■  CLTYPE 


1  • 

Cirrus 

5  - 

A1 tocuaulus 

2  - 

Cirrostratus 

6  • 

Cumulus 

3  - 

Stratus 

7  - 

St ratocumulus 

4 

A1 tost  ratus 

8  - 

Cumulonimbus 

Weather 

Status  ■  WESTAT 

A  - 

Hail 

K  * 

Smoke 

BN- 

Blowing  Sand 

L  - 

Drizzle 

D  - 

Dust 

RW- 

Rain  Showers 

F  - 

Fog 

S  - 

Snow 

GF- 

Ground  Fog 

T  - 

Thunderstorm 

H  - 

Haze 

ZL- 

Freezing  Drizzle 

Precipitation  -  PRECIP 

Numerical  -  Amount  of  precipitation 

in  tenths  of  an  inch. 

OR 


H  -  Heavy 
L  -  Light 
M  -  Moderate 
S  -  Slight 
T  -  Trace 

Ground  Visibility  -  GRDVIS 

1  -  Excellent 

2  -  Good 

3  -  Fair 

4  -  Poor 
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Moon  Phase  -  MUNPHA 


A. 9 

1  -  One-quarter  full 

2  •  One-half  full 

5  -  Three-quarter  full 

4  ■  Full 

5  -  No  moon 

A.  10  Foliage  -  FPL  1  AG 

1  -  Open 

2  -  Light 

3  -  Medium 

4  -  Dense 
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updating  as  more  dats  are  complied.  Contained  vlthln  Is  a  complete  program 
documentation  Including  step  by  step  operating  procedures  and  logic  flow  diagrams. 
Tn  addition,  the  data  reduction  methodology  and  program  test  results  are  described. 
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